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Abstract 


This  technical  report  analyses,  theoretically,  the  capability  of  ground  mov¬ 
ing  target  detection  based  on  the  SAR  along-track  interferometric  phase. 
The  probability  density  function  of  the  interferometric  phase  is  derived  for 
the  general  case  that  moving  target  signals  are  superimposed  upon  the 
stationary  clutter.  This  theoretical  analysis  is  an  extension  of  the  known 
statistics  for  the  clutter-only  case.  Several  target  models  are  proposed  that 
depend  on  one  hand  on  the  spatial  dimension  of  the  moving  target  com¬ 
pared  to  the  multilook  resolution  cell  size,  and  on  the  other  hand  on  the 
backscattering  type,  i.e.,  deterministic  or  random  target  signal.  The  derived 
density  functions  provide  the  means  to  quantify  the  performance  limits  (by 
determining  the  receiver  operating  characteristics)  for  any  set  of  system  and 
target  parameters,  such  as  false  alarm  rate,  signal-to-clutter  ratios  (SCR) 
and  target  velocities.  Hence,  deeper  insight  is  provided  into  results  ex¬ 
pected  from  upcoming  SAR/GMTI  systems  such  as  the  Radarsat2  Modex 
experiment  or  the  enhanced  airborne  Aurora  CP-140  SpotSAR  system. 

Resume 


Le  present  rapport  technique  expose  une  analyse  theorique  sur  la  capacite 
de  la  detection  des  cibles  mobiles  basee  sur  l’interferometrie  longitudinale 
par  radar  a  antenne  synthetique  (RAS).  La  fonction  de  densite  de  proba¬ 
bility  de  la  phase  de  1’interferogramme  est  derivee  pour  le  cas  general  o 
les  signaux  des  cibles  mobiles  se  superposent  a  des  fouillis  stationnaires. 
Cette  analyse  theorique  est  la  continuation  statistique  du  cas  o  il  n’y  a 
que  du  fouillis.  Plusieurs  modeles  sont  proposes  selon,  d’une  part,  la  di¬ 
mension  spatiale  de  la  cible  mobile  par  rapport  a  la  taille  de  la  cellule  de 
resolution  a  multiples  facettes  et,  d’autre  part,  le  type  de  retrodiffusion  (c.- 
a-d.,  signaux  deterministes  ou  signaux  de  cibles  aleatoires).  Les  fonctions 
de  densite  derivees  offrent  le  moyen  de  quantifier  les  limites  de  performance 
(en  determinant  les  caracteristiques  de  fonctionnement  du  recepteur)  pour 
tout  ensemble  de  parametres  de  systemes  et  de  cibles,  tels  que  le  taux 
de  fausse  alerte,  les  rapports  signal/fouillis  et  les  vitesses  anticipees  des 
cibles.  Ainsi,  on  comprendra  mieux  les  resultats  qu’on  attend  des  prochains 
systemes  RAS/GMTI,  notamment  de  l’experience  sur  le  Radarsat2  Modex 
ou  le  systeme  perfectionne  SpotSAR  a  bord  du  CP-140  Aurora. 


DRDC  Ottawa  TR  2002-084 


i 


DRDC  Ottawa  TR  2002-084 


Executive  summary 


In  space-borne  and  airborne  radar  surveillance  it  is  desirable  to  be  able 
to  detect  military  moving  targets  such  as  tanks,  wheeled  vehicles  or  ships 
within  strong  ground  clutter.  One  possible  way  to  achieve  this  is  with  Syn¬ 
thetic  Aperture  Radar  Along-Track  Interferometry  (SAR-ATI),  where  two 
channels  observe  the  same  scene  at  different  times.  For  stationary  terrain 
the  two  channel  signals  are  identical  and  can  be  canceled  out  (clutter  sup¬ 
pression)  by  computing  the  phase  difference,  i.e.,  the  interferogram,  leaving 
primarily  the  moving  targets  in  the  differential  data. 

In  the  past,  SAR-ATI  was  mainly  treated  in  the  context  of  measurement 
of  sea  surface  currents,  though  the  capability  for  sensitive  ground  moving 
target  indication  (GMTI)  is  obvious.  Although  the  effectiveness  of  ATI 
for  GMTI  has  been  recently  evaluated  with  measured  data  from  several 
flight  experiments  no  thorough  theoretical  discussion  has  been  performed 
so  far.  Precise  knowledge  of  the  interferogram’s  phase  statistics  is  crucial 
for  the  development  of  statistically  based  detector  tests  for  distinguishing 
the  moving  targets  from  the  clutter. 

The  purpose  of  this  technical  report  is  to  investigate  the  probability  density 
function  of  the  SAR  interferogram’s  phase  when  target  signals  are  superim¬ 
posed  upon  the  clutter  echoes.  Depending  on  the  spatial  dimension  of  the 
moving  target  compared  to  the  SAR  resolution  cell  size  and  the  backscat- 
tering  type  (i.e.,  deterministic  or  random  target  signal),  several  models  are 
proposed.  The  derived  density  functions  for  most  models  are  graphically 
illustrated  and  numerically  confirmed  with  simulated  interferometric  SAR 
data. 

This  analysis  provides  the  means  to  quantify  the  performance  limits  of  fu¬ 
ture  SAR  systems  for  arbitrary  sets  of  system  and  target  parameters,  such 
as  false  alarm  rates,  signal-to-clutter  ratios  (SCR)  and  anticipated  target 
velocities.  Deeper  insight  can  be  gained  into  results  expected  from  upcom¬ 
ing  SAR/GMTI  systems  particularly  the  Radarsat2  GMTI  experiment  and 
the  enhanced  airborne  Aurora  CP-140  SpotSAR  system. 

One  remarkable  result  for  a  typical  space-borne  system  is  that  a  target, 
with  radial  velocity  of  at  least  =  6  m/s  and  with  signal  power  at  least 
three  times  the  surrounding  background  (SCR  =  5  dB)  can  be  detected 
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with  better  than  98%  probability  whereas,  on  average,  only  one  false  alarm 
occurs  in  every  10000  SAR  image  pixels. 

Christoph  H.  Gierull.  2002.  Moving  Target  Detection  with  Along-Track  SAR 
Interferometry.  DRDC  Ottawa  TR  2002-084.  Defence  R&D  Canada  -  Ottawa. 
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Sommaire 


Dans  la  surveillance  au  moyen  de  radars  aerospatiaux  et  aeroportes,  il  est 
souhaitable  de  pouvoir  detecter  en  presence  de  forts  fouillis  de  sol  des  cibles 
mobiles  militaires,  comrne  des  chars,  des  vehicules  a  roues  ou  des  navires. 
Un  moyen  possible  de  le  faire  est  d’avoir  recours  a  l’interferometrie  longi- 
tudinale  par  radar  a  antenne  synthetique  (RAS),  qui  permet  d’observer  la 
meme  scene  dans  deux  canaux  a  des  moments  different^.  En  terrain  fixe,  les 
signaux  des  deux  canaux  sont  identiques  et  il  est  possible  de  les  annuler  (sup¬ 
pression  de  fouillis)  en  calculant  le  dephasage,  c’est-a-dire  l’interferogramme, 
ce  qui  laisse  surtout  les  cibles  mobiles  dans  les  donnees  different ielles. 

Par  le  passe,  l’interferometrie  longitudinale  par  RAS  a  ete  consideree  prin- 
cipalement  dans  le  contexte  de  la  mesure  des  courants  de  surface  en  mer, 
meme  s’il  y  a  manifestement  une  capacite  d’indication  de  cibles  terrestres 
mobiles  (GMTI)  precise.  Meme  si  on  a  recemment  evalue  l’efhcacite  de 
l’interferometrie  longitudinale  pour  la  GMTI  au  moyen  de  donnees  mesurees 
lors  de  plusieurs  experiences  en  vol,  il  n’y  a  pas  encore  eu  d’examen  theorique 
approfondi.  Il  est  crucial  de  connaitre  avec  precision  les  donnees  des  phases 
de  l’interferogramme  pour  developper  des  essais  de  detection  fondes  sur  des 
statistiques  dans  le  but  de  distinguer  entre  les  cibles  mobiles  et  le  fouillis. 

Le  present  rapport  technique  a  pour  but  d’etudier  la  fonction  de  densite  de 
probability  de  la  phase  de  l’interferogramme  du  RAS  lorsque  les  signaux  de 
cibles  sont  superposes  a  des  echos  de  fouillis.  Selon  les  dimensions  spatiales 
de  la  cible  mobile  par  rapport  a  la  taille  de  la  cellule  de  resolution  du 
RAS  et  au  type  de  retrodiffusion  (c.’est-a-dire  des  signaux  deterministes 
ou  des  signaux  de  cibles  aleatoires),  plusieurs  modeles  sont  proposes.  Les 
fonctions  de  densite  derivecs  pour  la  plupart  des  modeles  sont  illustrees 
graphiquement  et  confirmees  numeriquement  par  des  donnees  simulees  de 
RAS  interferometriques. 

La  presente  analyse  fournit  le  moyen  de  quantifier  les  limites  de  performance 
des  futurs  systemes  RAS  pour  des  ensembles  arbitrages  de  parametres  de 
systeme  et  de  cible,  comme  les  taux  de  fausse  alarme,  les  rapports  sig- 
nal/fouillis  et  les  vitesses  des  cibles  anticipees.  On  pourra  ainsi  rnieux  com- 
prendre  les  resultats  attendus  des  prochains  systemes  de  RAS/GMTI,  en 
particulier  ceux  qui  ressortiront.  de  1’experience  sur  le  dispositif  GMTI  de 
Radarsat  2  et  le  systeme  ameliore  de  detection  aeroporte  SpotSAR  a  bord 
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du  CP- 140  Aurora.  Un  resultat  remarquable  pour  un  systeme  aeroporte 
type  est  qu’il  est  possible  de  detecter  une  cible  ayant  une  vitesse  radiale 
d’au  moins  6  m/s  et  un  signal  d’une  puissance  d’au  moins  trois  fois  le  bruit 
ambiant  (rapport  signal/fouillis  de  5  dB)  avec  une  probability  superieure  a 
98  %,  alors  qu’en  moyenne,  il  n’y  a  qu’une  seule  fausse  alarme  a  tous  les  10 
000  pixels  d’ images  de  RAS. 

Christoph  H.  Gierull.  2002.  Detection  de  cibles  mobiles  au  moyen  de  l’interferometrie 
longitudinale  par  radar  a  antenne  synthetique  (RAS)  -  Analyse  theorique 
DRDC  Ottawa  TR  2002-084.  R&D  pour  la  defense  Canada  -  Ottawa. 
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1.  Introduction 


In  air-to-ground  radar  surveillance  it  is  desirable  to  be  able  to  detect  small 
moving  targets  within  strong  ground  clutter.  One  possible  way  to  achieve 
this  is  with  Along- Track  SAR  Interferometry  (ATI).  ATI  exploits  the  dif¬ 
ference  in  the  echoes  of  two  channels  observing  the  same  scene  at  different 
times.  These  channels  are  aligned  in  the  flight  direction  of  the  platform.  For 
stationary  terrain  the  two  channel  signals  are  identical  and  can  be  canceled 
out  (clutter  suppression)  by  computing  the  phase  difference,  i.e.,  the  inter- 
ferogram,  leaving  only  the  moving  targets  in  the  differential  data.  Precise 
knowledge  of  the  interferogram’s  phase  and  amplitude  statistics  is  crucial 
for  the  development  of  statistically  based  detector  tests  for  distinguishing 
the  moving  targets  from  clutter.  The  statistics  of  interferograms  for  dif¬ 
ferent  types  of  stationary  clutter  have  been  intensively  studied  in  the  past 
[1,  2,  3,  4], 

The  multilook  phase  difference  between  the  two  channel  outputs  Zi  is 

*  =  arg  ,  (1) 

where  n  denotes  the  number  of  independent  samples,  traditionally  called 
number  of  looks.  Under  the  assumption  of  Gaussian  distributed  clutter, 
the  density  function  was  shown  to  be 

=  r(n  +  1/2)  (1  —  p2)n  pcosip 

2v/7iT(n)  (1  —  p2  cos2  ip)n+l^2 

+  ---~2^  ~  2^i (n,  1;  1/2;  p 2  cos2  ip)  (2) 

for  — 7T  <  ip  <  7r  [1,  3].  Herein,  2Ai(-)  denotes  the  Gauss  hypergeometric 
function,  see  e  g.  [5].  The  pdf  of  eq.  (2)  depends  only  on  the  number  of 
looks  n  and  the  magnitude  of  the  complex  correlation  coefficient  pexp(d). 

In  the  past  SAR- ATI  has  been  mainly  treated  in  the  context  of  measurement 
of  sea  surface  currents  [6,  7],  though  the  capability  for  sensitive  ground 
moving  target  indication  (GMTI)  is  obvious.  Although  the  effectiveness 
of  ATI  for  GMTI  has  been  recently  evaluated  with  real  data  from  several 
experiments  [8,  9]  no  thoroughly  theoretic  discussion  has  yet  been  provided. 
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The  purpose  of  this  technical  report  is  to  investigate  the  probability  density 
function  of  the  SAR  interferogram  when  target  signals  are  superimposed 
upon  the  clutter  echoes.  Depending  on  the  spatial  scale  of  the  resolution 
cell,  this  density  function  is  used  to  define  detection  tests  for  slowly  mov¬ 
ing  targets.  In  a  previous  work  [4],  only  the  simplest  statistical  test  case, 
namely,  whether  the  contents  of  the  image  cell  in  question  were  target  only 
or  clutter  only,  was  considered.  The  underlying  assumption  was  that  the 
dimension  of  the  moving  targets  is  on  the  order  of  the  spatial  geometric 
resolution  of  the  SAR  and  that  the  backscatter  power  contained  in  the  cell 
only  resulted  from  the  desired  target,  i.e.,  the  clutter  (or  speckle)  power  was 
negligible.  Therefore,  only  the  density  function  of  the  clutter  was  required 
to  create  the  test,  i.e.,  to  determine  a  suitable  threshold  [4]. 

In  cases  where  the  resolution  cell  size  is  larger  than  the  target  size,  the 
clutter  power  can  no  longer  be  neglected.  Here,  the  hypothesis  “clutter 
only”  has  to  be  tested  against  the  alternative  “clutter  plus  target”.  This 
general  test  problem  represents  the  more  practical  case  that  a  cell  contains 
both  clutter  and  target  at  the  same  time.  In  this  context,  it  is  interest¬ 
ing  to  note  that  the  SAR  typical  effect  of  spatial  displacement  of  moving 
targets  within  the  SAR  image  substantiates  the  assumption  that  a  target 
signal  is  always  contaminated  with  clutter  power.  Even  with  the  terrain 
(clutter)  underneath  the  mover  masked  out  by  the  vehicle’s  body,  the  body 
is  imaged  displaced,  and  hence,  superimposed  upon  clutter  at  a  different 
“wrong”  location.  This  “wrong  clutter”  might  not  even  be  physically  rea¬ 
sonable,  e.g.,  the  moving  target  is  imaged  on  top  of  a  building  instead  of 
an  adjacent  street,  which  may  result  in  more  difficult  detection  due  to  a 
decreased  Signal-to-Clutter  ratio  (SCR). 

Depending  on  the  spatial  dimension  of  the  moving  target  compared  to  the 
multilook  cell  size  and  the  type  of  backscattering  model,  i.e.,  deterministic 
versus  random  target  echo,  several  models  are  theoretically  investigated  in 
the  subsequent  chapters.  The  derived  density  functions  for  most  models 
are  then  graphically  illustrated  and  numerically  confirmed  with  simulated 
interferometric  SAR  data.  Finally,  the  performance  limits  of  such  detectors 
are  determined  via  two  well  known  measures  of  quality,  i.e.,  the  probability 
of  detection  versus  varying  SCR  on  one  hand  and  versus  the  false  alarm 
rate  on  the  other  hand.  The  latter  plot  is  generally  known  as  the  receiver 
operating  characteristics  (ROC)  of  a  detector.  Both  measures  depend  also 
on  the  chosen  number  of  looks  and  on  the  resulting  interferometric  Doppler 
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phase  value,  i.e.,  the  radial  velocity  of  the  moving  target. 


1.1  Stationary  Clutter  Model 

Let  the  zero-mean  complex  signals  of  N  different  channels  define  the  column 
vector  Z_  —  [Z\,  -  ■  •  ,  Zn]T  ,  where  T  denotes  the  transpose  operator.  Under 
some  reasonable  assumptions  [4],  Z,  with 

Zi  —  (JiXi  V  i  =  l,---,N  and  X*  ~  A/4  (0, 1) ,  (3) 

can  be  modeled  as  a  multivariate  complex  Gaussian  random  vector  with 
density 

where  R  denotes  the  covariance  matrix. 


The  following  notational  convention  has  been  adopted  throughout  the  re¬ 
port:  Random  variables  are  capitalized,  their  realizations  are  represented 
by  the  corresponding  lower  case  letter.  Vectors  are  underlined  and  matrices 
are  written  bold.  The  superscript  H  denotes  complex  conjugate  transpose 
and  det(R)  the  determinant  of  R  [10].  In  the  case  of  SAR  interferometry 
usually  only  two  channels  are  involved  [11],  i.e.,  the  dimension  of  Z  is  cho¬ 
sen  to  be  N  =  2  as  shall  be  adopted  from  here  onwards  (without  loss  of 
generality),  i.e., 


R  =  eZZh  = 


Cn 

/CnC^e--70/ 


/CnC22^°  P 
C22 


The  off-diagonal  elements  Iii2  =  R21  —  EZj  Z|  of  the  covariance  matrix  R 
describe  the  complex  interferogram.  While  for  ac.ross-track  interferometry 
the  interferometric  phase  9  in  eq.  (4)  depends  on  the  illumination  geom¬ 
etry  and  the  terrain  elevation,  for  along-track  interferometry  this  value  is 
assumed  to  be  zero  ( 9  =  0),  i.e.,  the  clutter  is  assumed  to  be  stationary.  In 
cases  of  internal  clutter  motion,  such  as  sea  surface  currents,  9  is  a  function 
of  the  current  velocity. 


In  order  to  reduce  the  speckle,  polarimetric  and  interferometric  data  are  fre¬ 
quently  multilook  processed.  Multilook  interferometric  processing  requires 
averaging  several  independent  one-look  interferograms.  In  other  words  sev¬ 
eral  independent  one-look  sample  covariance  matrices  are  averaged.  The 
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n-look  sample  covariance  matrix  is  given  as 


(5) 

n  fc=i 

where  n  is  the  number  of  looks  and  zk  the  k- th  one-look  or  single-look 
sample  (also  called  snapshots).  The  random  matrix  A  =  nR  with 


Au  9*(Ai2)  +  j3?(Ai2) 

3^(Ai2)  —  i&(Ai2)  A22 

An  ae 

ae-M  A22 


(6) 


is  well-known  to  be  complex  Wishart-distributed  A  W2C  (ra,  R)  [10]  with 
probability  density 

u  (A)  =  exp  (~  tr{R~‘A})  ■  m 

for  n  =  2,3,3,...  Based  on  this  density  function,  the  multilook  phase- 
difference  (eq.  (2))  and  multilook  amplitude  density  functions  under  the 
hypothesis  “clutter  only”  have  been  derived  in  [1,  2,  3],  see  also  [4]. 


1.2  Inclusion  of  Channel  Errors  and  Imbalances 

In  practice,  interference  on  the  returned  echo  signals  from  different  sources 
is  unavoidable.  Such  interference  could  be,  for  instance,  thermal  noise  in 
the  channels,  non-ideal  receiver  hardware  or  uncompensated  range  migra¬ 
tion  in  the  SAR  processor.  The  implicit  study  of  the  impact  on  the  joint 
density  function  of  phase  and  magnitude  of  the  SAR  interferogram  has 
been  neglected  in  [1,  2,  3],  but  has  been  done  for  the  marginal  phase  den¬ 
sity  function  (2)  in  [12].  This  subsection  describes  how  to  introduce  such 
additional  terms  into  the  joint  interferogram’s  probability  density  function 
(pdf).  As  an  example  additional  thermal  noise  is  studied.  The  general  case 
with  arbitrary  transfer  functions  is  introduced  at  the  end  of  the  subsection. 

In  the  following,  an  additional  interference  TV  =  [TV!, . . . ,  Nn]t  is  modeled  as 
complex  normal  distributed  and  stochastically  independent  between  chan¬ 
nels  as  well  as  independent  from  the  clutter  and  the  moving  targets  but  with 
different  power  levels  of  for  i  =  1, . . . ,  TV,  i.e. ,  TV  Af$f  (0,  diag  (of, . . . , 
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a%)).  Under  this  assumptions,  the  covariance  matrix  in  eq.  (4)  has  to  be 
extended  to 


Kn  =  E  ZZH  +  E  NNh 

=  R+  diag((7i,...,(j^) 


which  replaces  R  in  eq.  (7)  for  N  =  2.  Since 

det  (Kjv)  =  (Cn  +  of)  ( C22  +  of)  —  p2CnC22 
=  Cj  1 C22  ((1  +  /ci)  (1  +  /C2)  —  p2)  , 


(8) 


where  K\  =  of/Cn  and  k2  =  o\jC22  denote  the  inverse  of  the  Clutter-to- 
Noise  power  Ratio  (CNR),  and 


tr  (K^'A)  = 


An  (1  +  /c2)  jC22  +  ^22  (1  +  /ci)  /Cn  —  2ap\J  C\\C22  cos  ip 


C11C22  [(1  +  /ci)  (1  +  /c2)  —  p2} 
the  density  function  of  A  becomes 

{AnA22  —  a2)71 


(9) 


/(A) 


7rr(n)r(n  -  1  )CftC&  [(1  +  m)  (1  +  /c2)  -  p2]" 

(1  +  K,2)An/Cn  +  (1  +  K\)A22/ C22  —  2ap  costp  /  y/CnC22 
(1  +  /Cl)  (1  +  /c2)  -  p2 

(AiiA22  —  a2) 


exp  — 


7rr(n)r(n-l)^1^(l-p2)n 

An/Cn  +  A22/C22  —  2apcosi!>/ \fc~nC 


•exp  -- 


^2 


1  -  P 


(10) 


for  A\\A22  >  a2.  Eq.  (10)  is  identical  in  form  to  the  stationary  clutter  in 
[2],  but  has  modified  variances  Cn  =  Cn  (1  +  /c*)  for  i  —  1,2  and  modified 

coherence  _ 

P  =  p/\/(l  +  /ci)  (1  +  /c2). 

Following  the  derivation  in  [2],  it  becomes  evident  that  the  joint  density 
function  of  the  interferogram’s  phase  and  magnitude  remains  unchanged  in 
form  and  that  any  additional  interference  only  alters  (reduces)  the  coherence 
between  the  channels. 


Generally,  different  transfer  functions  between  channels,  representing  differ¬ 
ent  non-ideal  data  acquisition  hardware  and  also  realistic  non-ideal  process¬ 
ing  system  software  must  be  considered.  According  to  the  system  theoretical 


DRDC  Ottawa  TR  2002-084 


5 


model  in  [12],  the  coherence  can  be  written  as  p  =  p/y/q i<?2  with 

P  =  J  Hn{u)Hl2(u)H*y(u)H*2(u)  dw 

and 

Qi  ~  j  \Hn(u)Hi2(uj)\‘2  du  +  Ki  J  |tfi2(u;)|2  dco 

where  Hn(u> )  represents  the  transfer  function  of  the  data  acquisition  and 
Haiuj)  denotes  the  transfer  function  of  the  SAR  processor  in  channel  i  G 
{1,2}.  Simple  analytical  models  for  the  different  transfer  functions  can  be 
used  to  investigate  the  effect  on  the  coherence  and  therewith  on  the  density 
function  [12]. 
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2.  Deterministic  Moving  Target  Model 


As  a  first  intuitive  model,  a  deterministic  moving  target  signal  s  =  j3\  1 
is  considered  as  superimposed  upon  the  stationary  clutter  returns  in  each 
channel.  Without  loss  of  generality,  the  Doppler  phase  of  targets  in  the  fore 
channel  is  set  to  zero.  The  real  amplitude  (3  determines  the  signal-to-clutter 
ratio  (SCR).  Since  n  single- look  cells  of  the  interferogram  are  summed,  the 
resulting  statistics  depend  on  both  the  number  of  moving  targets  (included 
in  the  entire  multilook  cell)  and  their  dimension  compared  to  the  multilook 
cell  size. 


2.1  Multiple  Movers  in  Multi-Look  Cell 


In  the  most  general  case,  each  single-look  cell  contains  a  moving  point  scat- 
terer  with  different  RCS’  and  different  radial  velocities.  This  scenario  might 
happen  when  observing  either  urban  areas  with  a  potentially  large  number 
of  movers  or  extended  targets  where  different  scattering  centres  may  have 
different  radial  velocity  components,  such  as  large  ships  rolling  and  pitching 
in  the  sea.  In  such  cases,  the  multilook  interferogram  is  given  as 

/  =  \  j>i(*) + #*))*  M*) + •  (u) 

k=  1 

Since  the  targets  are  assumed  deterministic,  the  model  (11)  is  equivalent  to 
I  —  1/n  Y2=\  zi(k)*Z2{k)  with  zi(k)  =  Z\{k)  +  f3(k)  and  z2(k)  =  z2(k)  + 
(3{k)  exp(j,d(k)).  The  composite  random  vectors  z(k)  =  [zi(k),  z2{k)Y  are 
mutual  independent  identical  complex  normal  distributed  with  expectation 
m(k )  =  0(k)[ l,eJ1?^]T  for  k  =  1, . . .  ,  n  and  covariance  matrix  R.  In  con¬ 
trast  to  eq.  (5),  the  random  matrix 

n 

A  =  nR  =  (12) 

fc=i 


is  now  noncentral  complex  Wishart-  distributed  A  ~  (n,  R,  ft)  with 

noncentrality  matrix 


n  =  r-[mmh, 

where 

r  m  p( 2)  •••  p(n) 

[  p{2)e^2)  . . . 


(13) 

(14) 
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Using  Annex  B,  the  density  function  of  A  is 


/( A)  = 


det  (A)w~2 

7rr(n)r(n  —  1)  det  (R)n  f 


- tr(fl) g—  tr ( R- 1  A)  oFl(n;fiR-1A), 


for  A  being  positive  definite,  which  is  in  abbreviated  form  written  as  A  >  0 
from  here  onwards.  As  this  density  function  contains  a  hypergeometric  func¬ 
tion  of  matrix  arguments,  analytical  analysis  is  exceedingly  complicated. 
However,  for  two  practical  special  cases  statistical  evaluation  becomes  pos¬ 
sible. 


2.2  Single  Mover  in  Multi-Look  Cell 

Regarding  high  resolution  SAR  systems  with  resolutions  in  the  meter-range, 
the  assumption  of  a  single  moving  target  within  a  multilook  cell  for  a  rea¬ 
sonable  number  of  looks  ~  10  seems  adequate.  This  mover  can  either  be 
a  single  point  scatterer  or  an  extended  target  consisting  of  several  scatter¬ 
ing  centres  with  approximatively  the  same  amplitudes  and  radial  velocities, 
such  as  a  tank  or  a  large  truck  etc.  The  underlying  test  problem  reduces  to 

n 

1  =  4-  0)*  (^2 (k)  +  0e’’d)  Target  phase  $  known  (15) 

fc= l 

0  e  IR 

Hypothesis  H  :  P  —  0 

Alternative  K  :  0  >  0. 

Under  the  assumption,  that  the  target  covers  only  q  <  n  single  look  cells, 
the  matrix  of  expectation  vectors  in  eq.  (13)  is  given  as 

M  =  /j[  j.]/, 

where  the  vector  q  has  q  components  equal  to  one  and  the  rest  zeros,  e.g., 

£=  [1,1,...,1,0,...,0]T 

In  this  case  the  matrix  product  MMff  (and  therewith  also  the  noncentrality 
parameter  f2  =  R“1MMJJ)  has  rank  one.  Denoting  the  signal  vector  as 
s  =  0[l,  eJl9]T,  the  noncentrality  matrix  can  be  written  as  =  qH~1ssH . 
Using  the  annotation  in  Annex  C.2,  the  density  function  of  A  can  be  reduced 
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to  involve  only  hypergeometric  functions  of  scalar  arguments.  Changing  the 
variable  B  in  eq.  (C.3)  to 

A  =  W_1BW“h  =  RiUBU^Rj, 


with  Jacobian  J  =  det(R)2,  det  (B)  =  det  (R)  1  det  (A)  and  tr(B)  = 
tr  (R-1A)  leads  to 


/(A) 


det(A)n  2  -tr(R-’A)  -A 

7rr(n)r(n  —  1)  det  (Rj 


0Fi{n]\K), 


(16) 


where 

k  =  ufR-^AR-5Ml 

and 

A  =  ufR-5MMHR-3u1. 

Further  it  is  easy  to  verify  that  the  first  eigenvector  ux  and  the  only  nonzero 
eigenvalue  A  of  the  matrix 


R-5MM"R-i  =  gR~5ss"R-5 


are 

R“5s 

u,  =  — ,  , 

and 

A  =  q^K^s. 


(IT) 

(18) 


Inserting  eqs.  (17)  and  (18)  into  eq.  (16)  finally  yields  the  density  function 
of  A 


/( A)  = 


det  (A)n  2  ,  r  n 

7rr(n)"r(n  — T)  det  (R)n  “‘’HdR"  A}) 

•  exp  (— ?sffR-1s)  oF\(n;  gsHR-1AR-1s). 


(19) 


The  equivalent  density  function  for  real  variables  (when  the  noncentrality 
matrix  has  rank  one)  was  originally  derived  by  Anderson  and  Girshick  [13, 
14]  where  they  applied  the  identity 

o*i(a;  X)  =  hVx)  (20) 

vx  v  ' 
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for  the  confluent  hypergeometric  function  [15].  /„(•)  denotes  the  modified 
Bessel  function  of  order  v.  After  inserting  eqs.  (4)  and  (6)  into  eq.  (19),  the 
argument  of  the  hypergeometric  function  yields 

B2 

gsHR-1AR-1s  =  q- — —  [(l  +  p2  -  2p cosd)  (An  +  A22) 

+  2 a  (cos  (1?  +  tj))  +  p2 cos (i9  —  ■0)  —  2pcos-0)j  .(21) 

In  order  to  get  the  interferogram’s  joint  density  function  of  phase  and  mag¬ 
nitude,  eq.  (19)  has  to  be  integrated  with  respect  to  An  and  A2 2-  The 
analytical  integration  is  extremely  difficult  and  leads  to  awkward  terms  for 
the  density;  therefore,  instead,  a  numerical  integration  has  been  performed 
in  Fig.  1  to  illustrate  the  joint  density  function.  The  pdf  has  been  computed 
with  the  number  of  looks  equal  to  two  and  a  coherence  of  p  =  0.95.  The 
target  phase  is  =  1.3  rad  and  the  SCR  has  been  chosen  to  0  dB. 

The  peak  of  the  pdf  is  seen  to  be  centered  between  the  target  and  clutter 
phases,  and  possesses  nonsymmetric  tails  in  both  magnitude  and  phase. 
The  migration  of  the  density  function  from  the  clutter  phase  at  zero  towards 
the  target  phase  for  increasing  SCR  is  demonstrated  in  Fig.  2,  where  the 
marginal  density  function  of  the  phase  is  plotted.  When  the  SCR  tends  to 
infinity,  the  phase  density  function  becomes  a  delta-function  at  the  clutter 
Doppler  phase  because  the  random  clutter  component  is  negligible. 
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Nlooks  =  2,  p  =  0.95,  8  exp(j  9)  =  1  exp(j  1 .3) 


0 


2. 


Phase  [rad] 


Figure  1  :  Joint  probability  density  function  of  the  SAR  interferogram  including  a 
deterministic  target  with  phase  1.3  rad  and  SCR  of  0  dB,  Number  of 
looks  n  =  2. 


Phase  [rad] 


Figure  2  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  deterministic  target  with  phase  1.3  rad  and  varying  SCR.  Number  of 
looks  n  =  2. 
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2.3  Large  Number  of  Looks 

A  complete  analytical  description  of  the  joint  density  function  is  feasible 
for  the  case  where  the  number  of  looks  is  large.  After  multiplying  out  the 
terms  in  eq.  (15) 

I=lYl  W*)**&(*)  +  ^ ^ Z  HkY&*  +  faW) . 

fc=i  fc= i 

it  can  be  seen  that  the  second  sum  (of  independent  random  variables  with 
zero  mean)  tends  to  zero  for  increasing  n.  Hence  the  test  problem  can  be 
approximated  to  a  superposition  of  the  target  signal  on  the  interferogram 
rather  than  on  the  separate  channel  outputs: 

1  n 

I  =  —  (zi(k)*Z2(k)  +  0e J!?)  Target  phase  ft  known  (22) 

k= 1 

0eR 

Hypothesis  H  :  0  =  0 

Alternative  K  :  0  >  0. 


For  notational  convenience,  the  parameter  02  describing  the  amplitude  of 
the  target  signal  has  been  replaced  by  0  without  loss  of  generality.  To 
calculate  the  interferogram’s  statistics  under  the  alternative  in  eq.  (22),  the 
derivation  of  Lee  et.  al.  [2]  for  the  stationary  clutter  has  to  be  extended. 
Let 


S  =  ssH 


0  0ej*  ' 
0e~j*  0 


be  the  deterministic  matrix  describing  the  moving  target  signal,  such  that 


B  =  A  +  S  =  nR  +  S  =  ^2zkz^  +  S.  (23) 

fc=i 


The  interesting  random  variables,  describing  the  interferogram  of  clutter 
plus  target,  are  again  the  off-diagonal  elements  of  B.  The  matrix  B  is 
positive  definite  for  n  >  2.  The  density  function  of  B  can  therefore  be 
calculated  as 


/b(B)=  /a(B-S) 

det(B  -  S)n"2 


7rT  (n)  T(n  —  1)  det  (R) 


n  exp  (—  tr{R  1  (B  —  S)})  .  (24) 
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Using  eq.  (6),  the  determinant  in  eq.  (24)  is  given  as 


det(B  -  S)  =  (Au  -  0)  {A22  -  0)  -  (ae*  -  pej*)  ( ae~ ^  -  pe^) 

=  (An  —  0)  (A22  —  P) 

— a 2  —  02  +  2a/?cos  (ip  —  18)  >  0.  (25) 


Using  the  inverse  of  the  covariance  matrix 

1 


R 


-l 


C11U22  (1  —  P 2) 
the  trace  is  given  as 


C22  ~P\/CuC22 

— Pv/C'nC'22  Cn 


tr{R_1  (B  -  S)}  = 


(An  ~  P)  / Cn  +  (A22  —  P)  / C22 

1  -  p2 

2pQ/y/C,11C,22cos(^)  +  2pPj y/ CnC22  cos('d) 

1  —  p2 


(26) 


Changing  the  random  variables  from  (An  —  P)/Cn,{A22  —  /?)/C722 ,3F^( Ai2) , 
S(Ai2)  to  i?i,  02,p  =  a/ \jCnC22  introduces  the  Jacobian 

det(J)  =  1  /(C?nCt2v) 


(see  Appendix  A).  Now,  denoting  the  Signal-to-Clutter  ratio  (SCR)  as 
<5  =  p/ \fC  11^22  and  inserting  eq.  (25)  and  eq.  (26)  into  the  pdf  eq.  (24) 
yields 


rj,  if))  = 


(0i02  —  if  +  2  rj8  cos  {ip  —  1?)  —  82)n  2 
7r(l  —  p2)"T(n)r(n  —  1)  77 

cxp  ^  By  +  02  -  2p(pcosQ/;)  -8cos{8))^  ^ 


Note,  that  unlike  the  case  with  no  target,  /(0i,02,p,^)  is  still  a  function 
of  Cn  and  C22  because  of  <5.  The  pdf  of  magnitude  and  phase  is  obtained  by 
integrating  eq.  (27)  over  B\  and  02.  The  integration  domain  is  constrained 
by  B\B2  —  rf  +  2rj8cos{ip  —  1?)  —  82  >  0  because  the  matrix  B  is  positive 
definite.  Integrating  first  with  respect  to  B2  and  introducing  a  =  5  cos {ip-fl) 
and  =  <5  sin(-0  —  $)  yields 


f{Burj,ip) 


p0"  exp 


B 1  |  2p{rj  cos('f/,  i  —  d  cos(ip) 
1 —p2 


7r(l  —  p2)"T(n)r(n  —  1) 
(p  -  o)2  +  b 2 


(28) 


POO 

I  (17— n)2+b2 
B 1 


02 


01 


exp 


02 


1-P2 


d02 
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Using  the  integration  identity  from  Gradshteyn  and  Ryzhik  [15]  (p.  343, 
3.382,  2.)  we  deduce 


p5"  2  f  (p  -  a)2  +  b2\ 

7T(1  -  p2)r(n)  exp  V  (1-P2)5i  ) 

■  exp  (-  ( By  +  2p[pcos(-0)  -  5cos(i?)])  /(I  -  p2))  .  (29) 


Next,  we  integrate  eq.  (29)  with  respect  to  By 


f{Bi,r),ip) 


7r(l  —  p2)T(n) 


exp 


poo 

■  /  Br2 

Jo 


exp  - 


2p  (p  cos(^)  —  5  cos(i?)) 

1  —  p2 

(p  -  a)2  +  b2 


(30) 


By 


1  —  p2 


exp  - 


(1  -  p2)By 


dBy. 


Applying  the  integration  identity  from  Gradshteyn  and  Ryzhik  [15]  (p.  363, 
3.471,  9.)  gives 


f( _  ,/a  _  2p  ((p  —  a)2  +  b2)"2  ( 2 p  (pcos(^)  -5cos(t?)) 

7r(l  —  p2)r(n)  PV  1-P2 


•K»-i 


^^/(p-Q^  +  b2' 

v  1_p2 


(31) 


where  AT,^-)  denotes  the  modified  Bessel  function  of  order  v.  Taking  into 
account  that  the  sample  covariance  matrix  is  normalized  by  the  number  of 
looks,  i.e. , 


B  =  —  A  +  S 

n 


and 


/b(B)  =  /A  (nB  -  nS)  •  n, 
the  bivariate  density  function  is  finally  given  as 

2  nn+1p  ((p  —  6  cos  (ip  —  i?))2  +  52  sin  (ip  —  i?)2) 


n  — 1 
2 


f(v,tP)  = 


■  exp 


7r(  1  —  p2)T(n) 
f  2np  (pcos(V’)  —  5005(1?)) 


V 


1  —  p2 


/  2n-y/ (p  —  5  cos(^  —  1?))2  +  <52  sin(i/>  —  1?)2  \ 

( - w - )  ■ 


(32) 


The  pdf  (32)  is  illustrated  in  Fig.  3  for  a  number  of  looks  equal  three  and 
a  coherence  of  p  =  0.95.  The  target  phase  is  1?  =  1.3  rad  and  the  SCR 
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has  been  chosen  as  0  dB.  The  peak  of  the  pdf  is  centered  between  the 
target  phase  and  the  clutter  phase  at  zero  and  shows  a  large  tail  along 
magnitude  and  phase  (compare  also  Fig.  5).  After  increasing  the  number 
of  looks  to  n  =  9,  the  peak  of  the  pdf  moves  away  from  the  true  target 
phase  value  but  the  entire  pdf  becomes  more  symmetric,  Fig.  4.  Similar 
qualitative  behavior  can  be  seen  for  a  larger  SCR  of  20  dB  in  Fig.  5  and  6. 
In  order  to  verify  the  derivation,  the  marginal  density  function  of  the  phase 
is  compared  to  simulated  data  for  varying  SCR,  Fig.  7  and  8.  The  analytical 
integration  of  eq.  (32)  with  respect  to  the  magnitude  p  is  extremely  difficult 
or  even  impossible.  Hence  numerical  integration  was  used  to  calculate  the 
theoretical  phase  pdf.  The  perfect  agreement  of  theory  and  simulation  is 
recognized.  Again,  the  phase  pdf  tends  towards  a  delta  function  at  the  true 
target  phase  when  the  SCR  approaches  infinity. 


Nlooks  =  3,  p  =  0.95,  5  exp(j  0)  =  1  exp(j  1 .3) 


Phase  [rad] 


Figure  3  :  Joint  probability  density  function  of  the  SAR  interferogram  including  a 
deterministic  target  with  phase  1.3  rad  and  SCR  of  0  dB.  Number  of 
looks  is  n  =  3. 
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Nlooks  =  9,  p  =  0.95,  5  exp(j  8)  =  1  expfl  1 .3) 


0.2  0.4  0.6  0.8  1  1.2  1.4 

Phase  [rad] 


Figure  4  :  Joint  probability  density  function  of  the  SAR  interferogram  including  a 
deterministic  target  with  phase  1.3  rad  and  SCR  of  0  dB.  Number  of 
looks  n  =  9. 


Nlooks  =  3,  p  =  0.95,  5  exp(j  8)  =  10  exp(j  1 .3) 


Phase  [rad] 


Figure  5  :  Joint  probability  density  function  of  the  SAR  interferogram  including 
a  deterministic  target  with  phase  1.3  rad  and  SCR  of  20  dB.  Number 
of  looks  n  =  3. 
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Normalized  Magnitude 


Nlooks  =  6,  p  =  0.95,  8  exp(j  0)  =  10  exp(j  1.3) 


Phase  [rad] 


Figure  6  :  Joint  probability  density  function  of  the  SAR  interferogram  including 
a  deterministic  target  with  phase  1.3  rad  and  SCR  of  20  dB.  Number 
of  looks  n  =  6. 


Nlooks  =  3  ,  p  =  0.95 


Figure  7  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  deterministic  target  with  phase  1.3  rad.  Number  of  looks  n  =  3. 
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N looks  =  6  ,  p  =  0.95 


Figure  8  :  Phase  probability  density  function  of  the  5AR  interferogram  including 
a  deterministic  target  with  phase  1.3  rad.  Number  of  looks  n  =  6. 
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3.  Gaussian  Moving  Target  Model 


In  the  previous  test  problems  it  was  assumed  that  the  single  mover  ampli¬ 
tude  was  constant  from  one-look  cell  to  one-look  cell.  In  practise  this  will 
not  be  the  case  because  varying  propagation  conditions,  system  instabili¬ 
ties,  and  primarily,  variation  of  the  target  RCS  caused  by  changing  aspect 
angles  are  reasons  for  amplitude  fluctuation.  If  a  target  consists  of  many 
single  reflectors,  then  the  backscattered  signal  is  composed  of  many  single 
contributions  with  quasi  random  relative  phases.  Accordingly,  it  can  be 
approximated  as  complex  normal  distributed.  Over  short  periods  of  time, 
e.g.,  corresponding  to  the  fast-time  or  range  direction  in  SAR,  the  aspect 
angle  changes  only  a  little  and  the  amplitude  can  be  assumed  to  be  con¬ 
stant.  This  model  is  usually  named  the  Swerling  I-case  [16],  In  the  flight 
direction  of  the  SAR,  i.e. ,  slow-time  or  azimuth  direction,  the  pulses  are  far 
enough  apart  that  the  amplitudes  of  the  single  returns  can  be  considered  as 
stochastically  independent.  This  so-called  Swerling  II-case  is  investigated  in 
the  subsequent  sections.  Depending  on  the  geometric  resolution  compared 
to  the  target  dimensions  two  Swerling  II  applications  have  to  be  differenti¬ 
ated. 


3.1  Mover  Dimension  on  the  Order  of  the  Multi-Look 
Cell  Size 


Analogous  to  the  test  problem  in  eq.  (11)  the  alternative,  “clutter  plus 
target”,  can  be  modeled  as 

1  " 

/  =  -  ^2(zi(k)  +  Si(k))*  (Z2(k)  +  S2(k))  (33) 

fc=  i 


where  the  deterministic  target  signals  have  been  replaced  by  the  random 
variables  Si(k)  for  i  —  1,2.  As  mentioned  above,  the  composite  random 
vector  A  =  [Si,^]7  is  assumed  to  be  complex  normal  distributed  A  ~ 

J\f^  (0,  Ks)  with  expectation  zero  and  covariance  matrix 


K  S  =  P 


1  psej» 

pse-j»  1 


where  ps  denotes  the  coherence  of  the  moving  target  signals.  It  is  impor¬ 
tant  to  note  that,  in  the  absence  of  temporal  decorrelation,  the  coherence 
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between  the  two  channel  outputs  will  be  exactly  identical  for  the  target  and 
the  clutter.  However,  even  with  temporal  decorrelation,  the  target  coher¬ 
ence  can  either  be  larger  or  smaller  than  the  surrounding  clutter  coherence. 
A  perfectly  steady  vehicle,  for  instance,  smoothly  moving  over  a  terrain  of 
vegetation  might  have  a  larger  coherence  than  the  vegetation.  The  contrary, 
where  the  vehicle  is  rolling  and  bouncing  over  terrain  which  is  otherwise  per¬ 
fectly  stationary  leads  to  a  lower  coherence  for  the  target  compared  to  the 
surrounding  terrain. 


Since  the  clutter  and  the  target  are  mutually  stochastically  independent, 
the  sum  of  them  will  still  be  complex  normal  distributed.  The  resulting 
sample  covariance  matrix 


A  =  nR  =  +  + 

fc=  1 

n 

=  (34) 

k=  1 


with  vk  =  zk  +  sk  is  therefore  complex  Wishart  distributed  A  W^(n,R 
+Ks)  with  n  degrees  of  freedom.  The  composite  covariance  matrix  is  given 
as 


k5  +  r  = 


1  +  0  P  +  0ps 

p  +  0pse~^  1  +  0 


-  (1  +  0) 


y/p*+(Ppj+2IJpp.  cos  tf  jatanf  ) 


y/p*+Ppj+2fipp.  cos0  -jatanf  ) 

1+0 


l+P 


which  has  the  same  form  as  eq.  (4).  Hence,  the  joint  density  function  of  the 
phase  and  magnitude 


2rin+lr)n  (  2nr)pcos(i{)  —  tpo) 

-=^r  exP 


7rP(n)  (1  —  p2) 


1 


where  Kn  is  the  modified  Bessel  function  of  order  n  [3, 
changed  compared  to  the  stationary  clutter  case,  except 
p  and  the  phase  ip  are  replaced  by 


K, 


71—1 


2nr) 

[  —  p2 


(35) 

1,  17]  remains  un- 
that  the  coherence 


P  = 


a/p2  +  02pl  +  2  0pps  cos  '0 

1+0 


(36) 
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and 


^  -  ijj0  =  ijj  -  arctan  f  j  (37) 

\P  +  (3ps  cos'd  J 

The  corresponding  marginal  phase  density  function  is  given  in  eq.  (2). 

The  coherence  of  the  target  and  the  clutter  are  the  same  in  the  following 
numerical  examples;  The  effect  of  different  values  is  investigated  later.  For 
increasing  SCR  values,  the  phase  density  becomes  wider  than  the  original 
clutter  pdf  because  the  factor  ^/l  +  f32  +  2(3 cos'd / {1+ (3)  is  always  less  than 
one  if  d  ^  0,  assuming  that  ps  =  p. 

At  the  same  time  the  peak  of  the  pdf  migrates  towards  the  target  phase  value 
d.  When  the  SCR  tends  to  infinity  (p  — »  oo),  p  tends  to  p  and  V’o  towards  d, 
i.e. ,  the  phase  density  function  in  that  case  is  identical  to  the  original  clutter 
pdf  except  for  the  Doppler  shift  towards  the  target  frequency.  This  behavior 
of  the  joint  density  function  is  illustrated  in  Figs.  9-12  for  different  SCR’s 
and  number  of  looks  n,  and  for  the  marginal  phase  densities  in  Figs.  13-14. 
The  perfect  agreement  of  theory  and  simulation  can  be  recognized  again. 


Figure  9  :  Joint  probability  density  function  of  the  SAR  interferogram  including  a 
Gaussian  target  with  phase  1.3  rad  and  SCR  of  0  dB.  Number  of  looks 


71=1. 
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Nlooks  =  1,  p  =  0.95,  p  exp(j  9)  =  10  exp(j  1.3) 

).5 
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’  -3  -2  -1  0  1  2  3 

Phase  [rad] 

Figure  10  :  Joint  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  with  phase  1.3  rad  and  SCR  of  20  dB.  Number  of 
looks  n  =  1. 


Nlooks  =  3,p  =  0.95,  p  exp(j  0)  =  1  exp(j  1 .3) 
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-3-2-1  0  1  2  3 

Phase  [rad] 

Figure  11  :  Joint  probability  density  function  of  the  SAR  interferogram  including  a 
Gaussian  target  with  phase  1.3  rad  and  SCR  of  0  dB.  Number  of  looks 
n  =  3. 
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N looks  =  3,  p  =  0.95,  p  exp(j  0)  =  10  exp(j  1 .3) 


-3-2-1  0  1  2  3 

Phase  [rad] 

Figure  12  :  Joint  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  with  phase  1.3  rad  and  SCR  of  20  dB.  Number  of 
looks  n  =  3. 


Nlooks  =  1  ,  p  =  0.95 


-1  -0.5  0  0.5  1  1.5  2 


Phase  [rad] 

Figure  13  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  with  phase  1.3  rad  and  varying  SCR.  Number  of 
looks  n  =  1. 
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N  looks  =  9  ,  p  =  0.95 


Figure  14  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  with  phase  1.3  rad  and  varying  SCR.  Number  of 
looks  n  =  9. 


In  order  to  study  the  effects  of  different  correlation  values  for  the  target  and 
the  clutter,  two  simulation  have  been  carried  out.  Figure  15  demonstrates 
the  case  of  larger  coherence  of  the  moving  target  with  ps  =  0.99  compared 
to  p  —  0.95  for  the  surrounding  clutter.  The  solid  curves  represent  the 
density  functions  for  the  case  that  both  target  and  clutter  possess  identical 
correlation  p  =  ps  —  0.95.  One  can  see  that  the  mismatch  becomes  more 
pronounced  for  targets  with  larger  SCR,  i.e. ,  the  resulting  density  functions 
are  narrower  compared  to  the  case  of  identical  coherence.  When  the  corre¬ 
lation  of  the  moving  target  decreases  relative  to  that  of  the  clutter,  Fig.  16, 
the  opposite  effect  occurs;  The  density  functions  widen  for  larger  SCR  com¬ 
pared  to  the  identical  coherence  case.  However,  both  figures  also  indicate 
that  the  overall  impact  of  different  correlation  values  is  not  expected  to  be 
drastical  in  practice,  particularly  because  the  values  are  not  assumed  to 
vary  too  much  within  the  aperture  time  (i.e.,  the  time  for  the  aft  channel 
to  travel  to  the  previous  position  of  the  fore  channel).  On  the  contrary,  the 
coherence  values  are  expected  to  be  close  in  most  practical  cases. 
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Nlooks  =  9  ,  p  =  0.95 


Phase  [rad] 

Figure  15  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  with  phase  1.3  rad,  coherence  ps  =  0.99  and  varying 
SCR.  Number  of  looks  n  =  9. 


Nlooks  =  9  ,  p  =  0.95 


Figure  16  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  with  phase  1.3  rad,  coherence  ps  =  0.91  and  varying 
SCR.  Number  of  looks  n  =  9. 
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3.2  Mover  Dimension  on  the  Order  of  the  Single- 
Look  Cell  Size 


In  cases  where  the  moving  target  dimension  is  smaller  than  the  multilook 
cell  size,  the  model  in  eq.  (33)  is  no  longer  adequate.  Instead,  we  get  for 
the  sample  covariance  matrix 

l  n 

A  =  tek  +  Sk)  te*  +  S-kf  +  zkz% 

fc=l  k=l+ 1 

=  B  +  C  for  l  e  {2,...,n}, 


when  it  is  assumed  that  the  mover  is  only  contained  in  the  first  l  one- 
look  cells.  Therefore,  the  random  matrix  A  consists  of  a  sum  of  two  inde¬ 
pendent  complex  Wishart  distributed  matrices  B  (/,R  +  K5)  and 

C  ~  (n  —  l,  R).  Assuming  that  there  exist  density  functions  for  B  and 
C,  i.e.,  3  <  l  <  n  —  1,  the  density  function  of  A  is  given  as  the  integral 


/(A) 


L 


/b(B)/c(A  —  B)  dB  (38) 

B>0,A— B>0 

det  (R)_(n_i+1)  det  (R  +  Ks)-(i-1)  .  ^  \ \ 

7T2r(n  -  l  +  1  )T(n  -  l)T{l  -  1  )T{1  -  2)  CXP  ^  tr  +  Ks)  ^ 


L 


det(B)n_/_1det(A-B)("3e“tr((R“1"(R+Ks)  1)B)  dB. 


B>0,A-B>0 


To  the  best  knowledge  of  the  author,  such  a  density  function  has  not  yet 
been  derived  in  closed  form. 


However,  to  investigate  the  difference  between  the  phase  densities  in  section 
3.1  and  3.2,  a  numerical  study  has  been  performed.  Figure  17  shows  the  re¬ 
sulting  phase  histograms  when  only  the  three  first  one-look  cells  contain  the 
moving  target.  Other  than  that,  all  other  parameters  are  chosen  identical 
to  those  used  in  Fig.  14.  Two  differences  can  be  recognized:  on  one  hand 
the  effective  SCR  is  only  a  third  of  the  original  SCR  (3/2>  (because  the  target 
covers  only  a  three  single-look  cells  out  of  nine),  and  on  the  other  hand  an 
increase  in  variance.  The  resulting  widening  of  the  density  functions  can  be 
approximatively  represented  by  a  loss  in  the  effective  number  of  looks.  In 
Fig.  17,  for  instance,  an  effective  number  of  looks  of  n  =  6  is  assumed  for 
the  superimposed  theoretical  density  functions.  While  this  loss  can  partly 
be  compensated  by  choosing  a  larger  number  of  looks  for  the  processing, 
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the  reduction  in  SCR  is  inherent  and  will  of  course  reduce  the  probability 
of  detection,  compare  also  section  4. 


Figure  17  :  Phase  probability  density  function  of  the  SAR  interferogram  including 
a  Gaussian  target  in  only  three  one-look  cells  with  phase  1.3  rad  and 
varying  SCR.  Number  of  looks  n  —  9.  Theoretical  curves  are  calculated 
with  (3/3  and  n  =  6. 
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4.  Measures  of  Quality 


28 


In  order  to  study  and  compare  the  performance  of  the  different  models  for 
the  detectors,  two  well  known  measures  of  quality,  i.e. ,  the  probability  of 
detection  Pd  versus  SCR,  and  Pj,  versus  the  false  alarm  rate  are  studied.  As 
a  worst  case  scenario,  only  the  results  for  the  Gaussian  random  model  are 
exemplarily  introduced,  but  they  could  easily  be  computed  for  the  deter¬ 
ministic  target  model  as  well.  The  coherence  of  the  target  and  clutter  are 
assumed  to  be  identical. 


The  probability  of  detection  is  defined  as 

Pd(0  =  I 

r  T(n  +  1/2)  (1  -  p2)n  pcos(ip  -  ip0) 

h  2v/(7r)r(n)  (1  -  p2cos2('0  - 'i/’o))n+1/2 

+  -^  J*  ^  2-^1  (n,  1;  1/2;  p 2  cos2(^  -  ip 0))  dip  (39) 

Z7T 

where 

p  =  p\J  1  +  /32  +  2,3  cos  3/(1  +  f3) 

and 

ipo  =  arctan(/3sin$/(l  4-/?  cost?)). 


The  threshold  £  must  be  determined  in  such  a  way  that  the  probability  of 
false  alarm  (or  false  alarm  rate) 

/7T 

Mfr,  o,o)  dip 

T(n  +  1/2)  (1  —  p2)n  pcosip 

2y/ln)r(n)  (1  —  p2  cos2  ip)n+lP 
/-.  2\n 

4 — — ^ - — iFi{n.  1;  1/2;  p2  cos2  ip)  dip  <e  (40) 
2ir 

is  always  less  or  equal  than  a  given  level  e,  i.e.,  £  =  Fpl  (1  —  e),  where 
Fq,(-)  denotes  the  probability  distribution  function  of  4'.  In  the  following 
the  integrals  in  eq.  (39)  and  eq.  (40)  are  numerically  evaluated  for  different 
sets  of  target  parameters. 
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4.1  Probability  of  Detection  for  varying  SCR 

In  Figure  18  the  dependence  of  the  probability  of  detection  Pd  on  the  SCR 
and  the  chosen  number  of  looks  is  investigated.  The  coherence  is  p  =  0.95 
and  the  target  Doppler  phase  is  held  constant  at  74°  (1.3  rad).  It  can  be 
recognized  that  moving  targets  cannot  be  detected  even  for  very  large  SCR 
when  the  number  of  looks  is  small  n  <  4.  The  reason  for  this  deficiency 
is  the  relatively  broad  phase  distribution  for  small  number  of  looks  so  that 
the  threshold  (according  to  the  chosen  value  e  =  le~4)  becomes  too  large, 
especially  larger  than  d  =  1.3  rad  (see  Table  1).  On  the  other  hand,  for 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

vArad\ 

3.123 

2.840 

1.450 

0.846 

0.664 

0.548 

0.448 

0.413 

0.379 

0.341 

Table  1:  Phase  thresholds  for  varying  number  of  looks 


larger  values  of  n  (n  >  7)  only  a  small  additional  gain  in  performance  can 
be  achieved. 


p(  =  o.oooi,  a  =  74  deg 


Figure  18  :  Probability  of  detection  versus  SCR  for  varying  number  of  looks. 


In  Fig.  19  it  is  seen  that  by  holding  the  number  of  looks  constant  at  n  = 
9,  and  varying  the  Doppler  phase  d  from  zero  to  90°  a  minimum  target 
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radial  velocity  corresponding  to  about  d  =  40°  is  necessary  in  order  to  be 
able  to  detect  the  targets  for  reasonable  SCR.  According  to  Table  1,  the 
phase  threshold  was  set  to  rj^  =  0.379  for  n  =  9.  Again,  after  exceeding  a 
certain  target  velocity,  no  significant  additional  gain  in  detection  probability 
is  achieved.  Assuming  that  the  two  receiving  antennas  are  separated  in 
azimuth  by  the  distance  d,  the  interferometric  phase  can  be  a  approximated 


by  [18] 


A^>  =  ^ 

A  Vn 


or 


a  jl  A  Va 

= 


where  A  represents  the  wavelength,  vr  the  radial  velocity  of  the  target  and 
va  the  velocity  of  the  antenna.  Inserting  for  example,  typical  space-borne 
parameters:  A  =  5.6  cm,  va  =  7.5  km/s  and  d  =  7.5  m,  yields  for  A 0min  = 
40°  a  corresponding  minimum  target  velocity  of  tyirnin  =  6.3  m/s  given  a 


SCR  of  5  dB. 


Inserting,  for  instance,  some  realistic  airborne  parameters  for  along-track 
interferometry,  such  as  va  =  100  m/s  and  d  =  0.5  m,  yields  an  even  more 
sensitive  system  with  urimin  =  1-23  m/s. 


P,  =  0.0001 ,  Nlooks  =  9 


Figure  19  :  Probability  of  detection  versus  SCR  for  varying  target  Doppler  phases 
(target  velocities). 
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4.2  Receiver  Operating  Characteristics 

The  other  important  measure  of  quality  is  the  receiver  operating  character¬ 
istic  (ROC),  i.e.,  Pd  as  a  function  of  Pf  for  fixed  SCR.  In  Fig.  20  the  ROC 
is  shown  for  the  considered  test.  The  number  of  looks  was  chosen  as  n  =  5 
and  the  target  Doppler  phase  as  d  =  1.3  rad.  For  vanishing  SCR,  Pd  =  Pf 
holds;  the  corresponding  curve  is  the  diagonal  through  the  ROC-diagram. 
For  increasing  SCR  the  curves  bend  upwards  under  retention  of  the  points 
(0,0)  and  (1, 1).  A  best  test  for  the  alternative  is  characterized  in  such  a 
way  that  no  other  test  exists  whose  ROC  exceeds  that  of  the  best  test  at 
any  position.  The  ROC  provides  insight  into  the  minimum  allowable  Pf  for 
a  certain  SCR  without  loosing  too  many  true  targets. 


Nlooks  =  5,  =  74  deg 


Figure  20  :  Probability  of  detection  versus  false  alarm  rate  for  varying  SCR. 
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5.  Summary  and  Conclusions 


In  this  report,  the  marginal  probability  density  function  of  the  SAR  interfer- 
ogram’s  phase  has  been  derived  analytically  in  the  case  where  moving  target 
signals  are  superimposed  upon  the  stationary  clutter.  Based  on  the  spatial 
dimension  of  the  moving  targets  compared  to  the  multi-look  cell  size,  sev¬ 
eral  models,  including  deterministic  as  well  as  Gaussian  distributed  target 
returns  have  been  investigated.  This  theoretical  analysis  is  an  extension  of 
the  known  statistics  for  the  clutter-only  case.  Most  theoretical  results  were 
demonstrated  and  confirmed  with  simulated  interferometric  phase  data. 

All  considered  models  have  in  common  that  the  peaks  of  their  corresponding 
phase  density  functions  migrate  from  zero  (stationary  clutter)  towards  the 
target  Doppler  phase  (introduced  by  its  radial  velocity)  for  increasing  SCR. 
One  major  difference  between  the  deterministic  and  random  model  is  that 
the  density  of  the  first  tends  towards  a  delta  function  when  the  SCR  — >  oo, 
whereas  the  pdf  of  the  latter  preserves  a  certain  width  (variance). 

In  order  to  study  and  compare  the  performance  of  the  different  models  in 
more  detail,  two  well  known  measures  of  quality,  i.e.,  the  probability  of 
detection  Pd  versus  varying  SCR  and  Pd  versus  the  false  alarm  rate  (re¬ 
ceiver  operating  characteristics)  were  studied  based  on  the  derived  density 
functions.  These  measures  provide  insight  into  expected  capabilities  from 
upcoming  SAR/GMTI  systems  particularly  Radarsat2  and  the  enhanced 
airborne  Aurora  CP-140  SpotSAR  system. 

One  main  conclusion  of  this  report  is  that  there  exists  a  quasi-optimum 
value  for  the  number  of  independent  looks.  A  good  choice  for  the  number 
of  looks  for  the  interferometric  processor  is  in  the  order  of  n  —  9.  Smaller 
numbers  than  n  =  4  prohibit  almost  any  detection,  larger  numbers  do  not 
imply  any  significant  gain.  In  contrast,  increasing  the  number  of  looks 
unnecessarily,  leads  to  a  loss  in  geometric  resolution. 

A  reasonable  choice  for  the  false  alarm  rate  could  be  Pf  =  le-4,  because  it 
would  correspond  to  a  single  false  alarm  per  100  x  100  SAR  interferogram 
pixels.  For  a  SAR  system  with  geometric  resolution  in  meter  range,  this 
would  roughly  lead  to  one  false  alarm  in  an  one  square- kilometer  (1km2). 

For  such  parameter  adjustments  one  would  expect  that  a  space-borne  sys- 
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tem  may  be  able  (greater  than  98  %  chance)  to  detect  targets  that  have 
a  SCR  of  >  5  dB  when  moving  faster  than  approximatively  20  km/h.  An 
airborne  SAR  system  would  be  about  five  times  more  sensitive,  i.e. ,  the 
minimum  radial  velocity  of  the  moving  targets  would  be  in  the  range  of  4 
km/h  for  the  same  system  parameters.  Higher  velocities  or  larger  SCR,  of 
course,  improve  the  probability  of  detection. 
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Annex 

Jacobian 
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The  Jacobian  associated  with  the  change  of  variables 
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Its  determinant  is 


or 


det (J) 
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Annex  B 

The  Noncentral  Complex  Wishart  Distri¬ 
bution 


The  noncentral  Wishart  distribution  generalizes  the  noncentral  x2  distri¬ 
bution  in  the  same  way  the  central  Wishart  distribution  generalizes  the  x2 
distribution;  and  the  noncentral  complex  Wishart  distribution  is  the  exten¬ 
sion  of  the  noncentral  Wishart  distribution  to  the  complex  plane.  The  real 
Wishart  distribution  is  discussed  in  detail  in  [19,  20];  the  complex  Wishart 
distribution  was  introduced  in  [10]  and  the  inverse  complex  Wishart  distri¬ 
bution  in  [21].  Although  the  complex  noncentral  case  that  we  require  is  not 
difficult  to  develop,  it  is  noticeably  absent  from  the  references  above. 

If  the  columns  of  the  complex  random  matrix  Z  6  C  ,  i.e.,  the  random 

77lX  1 

vectors  Z^k  €  CC  for  k  =  1, . . .  ,n  are  independent  complex  normal  dis¬ 
tributed  J\f^  {Ukk i  R) >  then  A  =  ZZH  is  said  to  have  a  complex  noncentral 
Wishart  distribution  with  n  degrees  of  freedom,  covariance  matrix  R,  and 
matrix  of  noncentrality  parameters 

Q  =  R'MM11, 

TTl  X  71  /T 

where  M  =  [m^ , . . . ,  mn\  €  C  .  We  will  write  that  A~V^(n,R,ft). 

When  n  <  m,  A  is  singular  and  the  A  ~  ( n ,  R,  fi)  distribution  does 

not  have  a  density  function.  For  n  >  m,  the  density  function  of  A  is  given 
as 

/a(A)=  rm(n)det(Rrdel(Armexp(-ttn) 

•  exp  (—  tr  (R_1  A))  0Ai(n;  S1R_1A)  A  >  0,  (B.l) 

where  fl  =  R^MM"  and  Tm(n)  -  r  (n  -  i  +  1).  0F\(a;  X) 

is  the  confluent  hypergeometric  function  of  matrix  argument,  which  involves 
series  of  zonal  polynomials.  These  functions  have  been  defined  for  real  vari¬ 
ables  in  [20]  and  were  extended  for  complex  variables  in  [22] . 
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Annex  C 

Generalized  Bartlett’s  Theorem  for  the 
Noncentral  Complex  Wishart  Distribution 


C.l  Standard  Normal  Distribution 

This  subsection  extents  the  generalized  Bartlett’s  decomposition  when  the 
noncentrality  matrix  has  rank  one  to  the  complex  case.  For  real  variables 
it  is  given  in  Theorem  10.3.8  of  [20]. 

mxn 

Let  the  columns  of  the  complex  random  matrix  Z  £  (L  ,  i.e.,  the  random 

771X  1 

vectors  Z_k  £  C  for  k  =  1, . . . ,  n  are  identical  complex  normal  distributed 

(2ZL  I)  with  m  =  [mi,  0, . . . ,  0}1 ,  so  that  A  =  ZZH  is  A  ~  (n,  I,  ft) 

with  ft  =  diag  (|mi|2, 0, . . . ,  0).  Further,  let  the  Cholesky  decomposition 
of  A  be  A  =  QQn,  where  the  lower  triangular  matrix  Q  has  real  positive 
elements  Qu  on  the  diagonal.  Then  the  elements  Qij  of  Q  are  all  statistically 
independent, 

•  \Q\\  is  noncentral  x2  distributed  with  2 n  degrees  of  freedom  and  non¬ 
centrality  parameter  |mi|2,  xL(lmi|2)> 

•  \Q\  for  i  t=  2, . . . ,  m  are  central  x2  distributed  with  2 (n  —  i  +  1)  degrees 
of  freedom,  X2(„-i+i)  and 

•  all  off-diagonal  elements  Qt]  for  i  /  j  are  standard  complex  normal 
distributed,^^  (0,1). 

Proof: 


With  fi  having  the  above  form,  the  density  of  A  in  eq.  (B.l)  is 
/a (A)  =  -  1  det(A)n~'7lexp(-|?n1|2)exp(-tr(A))  0Pi(n;  |mi|2^n). 

Since  A  =  QQ^  we  have 
.  tr  (A)  =  tr(QffQ)=E7<iW« 

•  det  (A)  =  n™i  Qi  and  An  =  Q2n. 
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The  Jacobian  is  given  as  J  =  Qii™  *+1^  [20])  so  that  the  joint  density 
function  of  the  Qtj  for  1  <  i  <  j  <  m  can  be  written  as 

m 

/(Q)  =  n^exP  (-««««) 

3<i 

(~Qu)  exp  (-|mi|2)  0Fi(n;  Imil2^) 

m 

which  is  the  product  of  the  marginal  density  function  for  the  elements  of  Q 
stated  in  the  theorem. 


C.2  Extension  to  General  Case 

Using  the  generalized  eigenvalue  decomposition  for  the  hermitian-definite 
problem  described,  e.g.,  in  [23]  any  arbitrary  complex  random  matrix  Z 

with  independent  A/jf  {mi)  R)  distributed  column  vectors  {%  =  1, . . . ,  n) 
can  be  transformed  into  the  form  stated  in  the  theorem  C.l  as  long  as  its 
matrix  of  expectation  vectors  M  =  [m-,  , . . .  ,mn]  has  rank  one.  If 

W  =  R"5U, 

where  U  are  the  eigenvectors  of  i.e. , 

R^MMffR'2  =  UAUH  =  U diag  (A,  0, ... ,  0)  XJH ,  (C.2) 


(because  MMW  has  rank  one),  then  the  random  matrix 

X  =  WHZ 


has  the  expectation 

with 


EX  =  WhM 


EXEXh  =  UffR-*MMffR“*U  =  A, 
which  follows  from  eq.  (C.2).  Now  let 

B  =  XXH  =  WHZZHW  =  WhAW, 
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then 


EB  =  EZZHW  =  W;/RW 

=  U//R-5RR§U  =  UUh  =  I. 

Therefore,  matrix  B  is  (n,  I,  diag  (A,  0, ,  0))  distributed  with  density 

1  ,  /■ ry\7l—m 

/(B)  =  --e  —  exp(-tr  (B))  exp  (—A)  0*i(ra;  XBn)  (C.3) 

where 

J5n  =  ufR“5BR“2«1  (C.4) 

and 

A  =  nfR-5MMHR-5Ml.  (C.5) 

The  vector  ux  denotes  the  first  eigenvector,  i.e.,  the  first  column  of  U. 
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